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Abstract. Border bases arise as a canonical generalization of Grobner bases. Unfortunately, so far only 
border bases containing a reduced Grobner basis can be computed. In this note we extend our previous 
work [4] to arbitrary order ideals and devise a polyhedral characterization of order ideals and border 
bases: order ideals that support a border basis correspond one-to-one to integral points of the order ideal 

■ polytope. In particular, we establish a crucial connection between the ideal and its combinatorial structure. 
O^l ' Based on this characterization we adapt the classical border basis algorithm to allow for computing border 

, bases for arbitrary order ideals, which are independent of term orderings. We also show that finding a 

^ ' preference-optimal order ideal that supports a border basis is NP-hard. 

oo . 

^^j- 1. Introduction 

<^ I In this note we extend our previous work [4] to arbitrary order ideals. For the sake of completeness 

and self-containment we include the necessary preliminaries from f4l and establish the connection 
whenever appropriate. Grobner bases have been essential in commutative algebra to make important 
operations on ideals such as intersection, membership test, elimination, projection, and many more 
computationally accessible allowing to perform actual computations with ideals. Recently it became 
clear though that Grobner bases are not always the best choice in order to perform such operations. In 
I particular when the actual ideals under consideration are inferred from measured data. 

^ ■ Although every Grobner basis with respect to a degree-compatible term ordering can be extended to 

Q I a border basis (see [12, p. 281ff]), not every border basis is an extension of a Grobner basis and thus 

•/^ ■ the former constitutes a broader class of bases. Even if an order ideal has the right cardinality, it does 

not necessarily support a border basis. Both cases were exemplified in | |12[ Example 6]. The border 
I basis algorithm as formulated in [|12 || , which in turn is a specification of Mourrain's generic algorithm 

II18II . computes border bases of zero-dimensional ideals. Unfortunately, it does so for order ideals 
^ ', supported by a degree-compatible term ordering only and the computed border basis is supported by a 

reduced Grobner basis. It therefore falls short to compute border bases for more general order ideals. 
In [p^2, Proposition 5], an alternative algorithm was presented which can potentially compute more 
rS I general border bases but requires a priori knowledge of the order ideal that might support the border 

■ basis. Thus the algorithm allows for computing a border bases once the order ideal is known but does 
not solve the characterization problem itself Further, the actual algorithm is heavily based on Grobner 
basis computations which is unsatisfactory as pointed out in [12 , p. 284]. 

In many cases though it can be advantageous to compute border bases for order ideals that do not 
necessarily stem from a degree-compatible term ordering. It has been an open question to character- 
ize the admissible order ideals of a given zero-dimensional ideal. This question has been answered for 
degree-compatible order ideals in [4] (not to confuse with degree-compatible term orderings!). We will 
extend our previous results to provide a full characterization of all admissible order ideals. The concept 
of border bases and in particular the border basis algorithm is rather old and (re-) appeared in differ- 
ent fields of mathematics including computer algebra, discrete optimization, logic, and cryptography 



2000 Mathematics Subject Classification. Primary: 13P10; 90C57; secondary: 65H10;12Y05; 90C27; 68R05. 
Key words and phrases, order ideal polytope, border bases, Grobner bases, combinatorial optimization. 
Research partially supported by Hungarian Scientific Research Fund, grant No. K 67928 and by the German Research 
Foundation (DFG) funded SFB 805. 

1 



independently under different names (see e.g.,[l3l [T7lll8in ). A coherent framework was introduced in 
EH nulls] and [16, Section 6.4]. 

In discrete mathematics and combinatorial optimization, Alon's polynomial method (cf [HKll]) and 
variants of the border basis algorithm have been used to formulate and test satisfiability of large classes 
of combinatorial problems (see e.g., [9, Section 2.3], [8, 7], and references contained therein). An- 
other link between border bases and the Sherali-Adams closure, a method for convexification used in 
combinatorial optimization, was recently investigated in II20II . Border bases also have applications in 
cryptography where sparse quadratic systems have to be solved. These systems arise when rewriting 
the S-boxes of crypto systems as polynomial equations. For example, the celebrated XL, XSL, Mutan- 
tXL attacks (see e.g., [14,i5j) are essentially equivalent to the reformulation-linearization-technique of 
Sherali and Adams. 

Our contribution. Due to its dependence on degree-compatible orderings, the classical computation 
of border bases can only provide border bases supported by a reduced Grobner basis. These border 
bases do represent only a rather small fraction of all possible border bases (see [4, Example 2.6]) and 
border bases that have a particularly nice vector space basis often cannot be obtained by a degree- 
compatible term ordering. We will free the computation of a border basis from the term ordering by 
applying the theory of combinatorial optimization and polyhedral combinatorics to the combinatorial 
structure of zero-dimensional ideals. Our contribution is the following: 

Polyhedral characterization of all border bases. We revisit the polyhedral characterization of border 
bases for degree-compatible order ideals introduced in []4[] and extend it to arbitrary border bases. As a 
result, we provide a complete, polyhedral characterization of all border bases of any zero-dimensional 
ideal I. We associate an order ideal polytope P to I whose integral points are in one-to-one correspon- 
dence with order ideals supporting a border basis of I. This explicitly establishes the link between the 
combinatorial structure of the basis of the factor space and the structure of the ideal. Whether an order 
ideal supports a border basis is solely determined by the combinatorial structure of the order ideal 
polytope. A related result for Grobner bases of the vanishing ideal of generic points was established in 

El- 

Computing arbitrary border bases. We extend the algorithm in [4] which in turn is based on the classical 
border basis algorithm (cf [12]) and show how to compute border bases for arbitrary order ideals using 
the order ideal polytope without relying on a specific degree-compatible term ordering. 
Computing optimal order ideals. We will show that computing an optimal order ideal supporting a 
border basis with respect to a linear preference function on the monomials is NP-hard in general. As 
we merely ask for a basis transformation this is surprising in some sense. Interestingly, the NP-hardness 



does not stem from the hardness of computing the L-stabilized span as the problem remains NP-hard, 



even in cases where the L-stabilized span is small enough to be determined efficiently as shown in our 
reduction. 

Computational feasibility. The computational feasibility of our method has been shown in ^4^ for 
degree-compatible order ideals and translates directly to the more general setting here (we have one 
rather large block rather than several small ones). 

Outline. We start with the necessary preliminaries in [Section 2[ In [Section 3 1 we introduce the order 
ideal polytope and establish a bijection between its integral points and border bases (for a given ideal 
7). In [Section 41 we show how to compute border bases for generic order ideals. Finally, in [Section 51 
we establish the NP-hardness of computing optimal order ideals. 

2. Preliminaries 

We consider a polynomial ring K[X] over the field K with indeterminates X = {x^, ...,x„}. For 
convenience we define x"^ := rije[n] -^j ' for m ^ N" and let T" :— {x^ \ m € N"} be the monoid 
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of terms. For any d € N we let T"^ := {x"^ e T" | ||m||i < d} be the set of monomials of total 

degree at most d. For a polynomial p e K[X] with p — ^l^^cii^'"' we define the support of p to be 
supp(p) := {x""' I i e [Z]} and similarly, for a set of polynomials P c K[X] we define the support of P 
to be supp(P) := IJ^gp supp(p). Given a term ordering a, the leading term LTj^(p) of a polynomial p is 
LTo-(p) t with t G supp(p) such that for all t' e supp(P) \ {t} we have t >o. t'; the leading coefficient 
LCo-(p) of p is the coefficient of YI^{p). We drop the index cr if the ordering is clear from the context. 
Recall that the degree of a polynomial p e i<"[X] is deg(p) := max^rngg^pp^p-) ||m||i. The leading form LF(p) 

of a polynomial p = Xi!=i 'J;^"'' ^ K\%\ is defined to be LF(p) = Xi!=i ||m ||=d '^i^'"' where d = deg(p), 
i.e., we single out the part with maximum degree. Both LF and LT generalize to sets in the obvious way, 
i.e., for a set of polynomials P we define LF(P) := {LF(p) | p e P} and LT(P) := {LT(p) | p € P}. 

In the following we will frequently switch between considering polynomials M = {pi,...,p^}, the 
generated ideal, and the generated vector space whose coordinates are indexed by the monomials in 
the support of M. We denote the ideal generated by M as (M)jf [x] ^^d the vector space generated by 
M as (M)jf. For n e N we define [n] := {l,...,n}. All other notation is standard as to be found in 
flU EJ; we have chosen the border basis specific notation to be similar to the one in [|12 | 1 where the 
border basis algorithm was first introduced in its current form. 

Central to our discussion will be the notion of an order ideal: 

Definition 2.1. Let ^ be a finite subset of T". If for all t € ^ and t' e T" such that | t we have 
t' € 0, i.e., is closed under factors, then we call an order ideal. Furthermore, the border d oi a 
non-empty order ideal is the set of terms d := {Xjt \ j e [n], t ^ 0} \ 0. As an exception, we set 
50 {1} for the empty order ideal. 

Recall that an ideal / c K[X] is zero-dimensional, if and only if K[1L]/I is finite dimensional. The 
^-border basis of a zero-dimensional ideal / is a special set of polynomials: 

Definition 2.2. Let = {ti,...,t^} be an order ideal with border d0 = {bi,...,by}. Further let 
^ = {gi, . . . , gy} Q K[X] be a (finite) set of polynomials and let / c K[X] be a zero-dimensional ideal. 
Then the set ^ is an 0-border basis of / if: 

(1) the polynomials in have the form gj = bj — ^f^;^ (^ij^i for j € [v] and a^j € K; 

(2) mK[x] = I; 

(3) JC[X] = 7 ® (0)k as vector spaces. 

If there exists an ^-border basis of / then the order ideal supports a border bcLsis of I. 

Note that the condition (^)j<:[x] = J is a consequence of ^ c /, the particular form of the elements in 
^, and K[X] = I ® {0)k- More precisely, (^#)k[x] + {^)k is closed under multiplication by the Xj, and 
hence it is an ideal. As it contains 1, it must be the whole ring, and hence = ^ by the modular 

law. See [|T3l Proposition 4.4.2] for another proof. 

In particular, an order ideal supports an ^-border basis of / if and only if i<C[X] = 7 ® {0)k- 
Moreover, for any given order ideal and ideal / the ^-border basis of / is unique as bj has a unique 
representation in K\X\ = I ® (0)fr for all j € [v]. Furthermore, as K[X] = I ® (0)fr it follows that 
\0\ = dim(^)jf is invariant for all choices of 0. The requirement for / being zero-dimensional is a 
consequence of this condition as well, as 5^ (and in consequence 0) as part of the output of our 
computation should be finite. Clearly, as a vector space, I has a degree filtration, i.e., I = UjgpjJ"' 
where := {p e / | deg(p) < i}. For a set of monomials we define 0^'^ := {m ^ \ deg(m) = i} 
(and similarly for < instead of =). Recall that an order ideal is called degree-compatible (to 7) if 

dim(^=% = dim(T^)^-dim^ 
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for all i e N. Note that not all order ideals supporting a border basis of a given zero-dimensional ideal 
/ have to be degree-compatible (see e.g., [4] for an example). Moreover, not every order ideal stems 
from a term ordering as shown in the follow^ing example: 

Example 2.3 (A general order ideal). The homogeneous ideal having a Grobner basis in the degree 
lexicographic term ordering 

2 2 2 4 

^2 "I" ■^1-^2 -"-1' -^1^2' "^2 

has the associated order ideal 

{l,Xi,X2,XiX2,xl,xl}. 

Via a basis change, w^e obtain another order ideal 

w^hich cannot come from a term ordering. The reason is that every Grobner basis of the ideal must 
contain a polynomial w^ith degree-tw^o leading term, which therefore must be a multiple of the first 
generator X2 + X1X2 + xf modulo higher-degree terms. The degree-two leading term cannot be X1X2 
because if e.g., x^ < X2 then the term X2 is larger then XiX2. So the leading term must be either x^ or 
X2, which excludes the order ideal. 

In fact, it turns out that border bases with order ideals coming from a (degree-compatible) term 
ordering do cover only a small fraction of all possible border bases as exemplified in [4, Example 2.6]. 

2.1. Computing stable spans. Without proofs, we recall the underlying stable span computation of 
the classical border basis algorithm introduced in [12J as it will serve as a basis for our algorithm. 
The interested reader is referred to [flSl lllll for a general introduction to border bases and to [12J 
in particular for an introduction of the border basis algorithm. The classical border basis algorithm 
calculates border bases of zero-dimensional ideals with respect to an order ideal & which is induced by 
a degree-compatible term ordering u by successively generating a vector space approximation of the 
ideal. These approximations are generated via the following vector space neighborhood extensions: 

Definition 2.4. (cf. [12^, Definition 7.1 and paragraph preceding Proposition 13]) Let V be a vector 
space. We define the neighborhood extension of V to be 

V+ --v + Vxi + --- + Vx^. 

For a finite set W of polynomials, its neighborhood extension is 

W+ = WUWxiU---UWx„. 

Note that for a given set of polynomials W such that {W)jf = V we have (W^)^ = = as 

multiplication with X; is a -vector space homomorphism. It thus suffices to perform the neighborhood 
extension on a set of generators WofV. 

Let F be a finite set of polynomials and let L c T" be an order ideal, then F n (L)jf = {/ e F | 
supp(/) c L}, i.e., F n (L)jf contains only those polynomials that lie in the vector space generated by L. 
Clearly, for L = T"^ we have n (L)jf = {F}^'^. Using the neighborhood extension we can define: 

Definition 2.5. (cf |]T2l Definition 10]) Let L be an order ideal and let F be a finite set of poljTiomials 
such that supp(F) c L. The set F is L-stabilized if (F+)j^ n {L)j^ - (F)^^. The L-stable span F^ of F is 
the smallest vector space G containing F satisfying G"*" n (L)jf = G. 

A straightforward construction of the L-stable span of F is to inductively define the following in- 
creasing sequence of vector spaces: 

Fo:-(F)j, and F^+i :- F+ n (L)j, for ?c> 0. 

The union [J]^>qFi^ is the L-stable span F^ of F. 
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The set L represents our computational universe and we will be in particular concerned with finite 
sets L c T". Note that L-stability is a property of the vector space and does not depend on the basis: 

Remark 2.6. The L-stable span of a finite set F depends only on the generated vector space (F)jf, as 

In the following we will explain how the L-stable span can be computed explicitly for L = T"^. We 
will use a modified version of Gaussian elimination as a tool, which allows us to extend a given basis V 
with a set W as described in the following: 

Lemma 2.7. 11121 Lemma 12] Let V — {v^, . . . , v^} Q K[X] \ {0} be a finite set of polynomials such that 
LT(v;) LT(vj) whenever i,j e [r] with i j and LC(v;) = 1 for all i e [r]. Further let G = {gi, ...,gs} 
be a finite set of polynomials. Then Algorithm 12.81 computes a finite set of polynomials W Q K [X] with 
LC(w) = 1 for all w e W, LT(ui) LT(u2) for any distinct Uj, U2^VUW, and {V U W)^ = {V U G)^. 
(V, W may be empty.) 

Algorithm 2.8 (Gaussian Elimination for polynomials — GaussEl). 

Input: V, G as in Lemma \2y\ 
Output: W QK[X] as in Lemma\Z7\ 

(1) Let H:=G and rj := 0. 

(2) IfH = <d then return W := {v^+i, . . . , v^+^} and stop. 

(3) Choose f and remove it from H. Let i := 1. 

(4) /// = Oori>r + ri then go to step 

(5) 7/LT(/) = LT(Vi) then replace f with f - LC(/) • V;. Set i := 1 and go to step @. 

(6) Set i:=i + l. Go to step @. 

(7) If f then put -q := -q + 1 and let v^^^ :=//LC(/). Go to step 1^. 

We can now compute the L-stable span using the Gaussian elimination algorithm l2.8t 

Lemma 2.9. |]T2l Proposition 13] Let L — T"^ and F QKlX] be a finite set of polynomials supported on 
L. Then Alsorithm \2A0\ computes a vector space basis V of F^ withpairwise different leading terms. 

Algorithm 2.10 (L-stable span computation — LStabSpan). 

Input: F, L as in Lemma \2.9\ 
Output: V as in Lemma \2.9\ 

(1) V :=GaussEl(0,f). 

(2) W' — GaussEl(V, V+ \ V). 

(3) W := {w e I supp(w) QL} = {w sW' \ deg(w) <d}. 

(4) If\W\ > setV :=¥ UW and go to step ©• 

(5) Return V. 

The rationale for computing a stable span approximation is due to the following proposition that 
serves as a criterion for testing whether an order ideal supports a border basis. 

Proposition 2.11. ['12', Proposition 16] Let L be an order ideal. Further let I be an L-stabilized generating 
vector subspace of a zero-dimensional ideal I c iC[X], i.e., I"*" n (L)jf = I and {l) i^^^^ =1- If is an order 
ideal such that (L)jf = 7 ® (^)jf and d& Q L then 6 supports a border basis of I. 

We obtain the following corollary: 
Corollary 2.12. Let I be an T^-^-stabilized vector space satisfying I + (t'^^_^)j^ = ('^"d)^- 
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Proof. We apply Proposition [2JT] with the choice L := T^,^, I := {l) j^^^^ and & := T^^ \LT(/) where the 
leading terms are with respect to any degree-compatible term ordering. Clearly, (t"^^^ = I ® {G)k- 
The condition / + = (^'"jj)^^ ensures that consists of monomials of degree less than 

d, so 3^ c T"^. Hence the proposition applies, and we obtain iC[X] = / © {6)j(^. Together with 

For a worst-case bound on d, we will use the dimension of i<C[X]/7. The necessary technical back- 
ground is the following lemma. 

Lemma 2.13. Let I he a zero-dimensional ideal o/iC[X], and let d := dimiC[X]/7. Then 

(1) 7^'i/7^d-l ^ 

(2) (T^d.i)^//-"-' =if[X]/7. 

Proo/ Choose a degree-compatible term ordering. The associated order ideal (as every order ideal of 
size d) contains monomials of degree less than d. This proves I + (t!1^_^)^ = iC[X], from which the 
statements easily follow via the modular law. □ 

The border basis algorithm decomposes into two components. The first one is the calculation of the 
L-stable span for L = T" ^ for sufficiently large d (this is the main 'work') and the second component 
is the extraction of a border basis via the final reduction algorithm. Effectively, after using any degree- 
compatible term ordering in order to compute the L-stable span one can choose a different ordering 
with respect to which the basis can be transformed. This approach is very much in spirit of the FGLM 
algorithm, Steinitz's exchange lemma, and related basis transformation procedures. In Remark 12.61 it 
is shown that L -stability does not depend on the basis and thus any sensible basis transformation that 
results in an admissible order ideal & is allowed. In the following we will characterize all admissible 
order ideals. 



3. The order ideal polytope 

We will now introduce the order ideal polytope P (a 0/ 1 polytope) that characterizes all admissible 
order ideals that support a border basis (for a given zero-dimensional ideal 7) in an abstract fashion 
independent of particular vector space bases for the stable span approximation. Its role will be crucial 



for the later computation (Algorithm 4.3) of border bases for general order ideals. We will first focus on 
its properties and structure, then in lSection 41 we will consider the computational aspect. We will show 
that the integral points of the order ideal polytope are in bijection with those order ideals that support 
a border basis of a given zero-dimensional ideal I (Theorem 13. 2D . In order to do so, we approach the 
problem from a polyhedral point of view to capture the intrinsic combinatorics of having to ensure that 
K[K] = I ® {0)k as vector spaces on the one hand and being an order ideal on the other hand. 
The role of the polyhedral description becomes prominent in lSection 41 when the directness of the sum 
/ © {0)k is rephrased in the language of matrices and vector space bases. 

Definition 3.1. Let 7 be a zero-dimensional ideal. Its order ideal polytope P{I) is given by the system of 



inequalities in Figure 3.1 



To obtain a finite dimensional polytope, we bound the degree of the monomials by dimK"[X]/J. This 
bound is large enough to contain the worst case. 

We are ready to relate the order ideal polytope with order ideals. From now on, let A(7) denote the 
set of order ideals (that support a border basis) of a zero-dimensional ideal 1. li s A(7) we call 
admissihle (for 7). 
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Figure 3.1. Order ideal polytope P(I) with d := dim 7C[X]// 

Theorem 3.2. Let I be a zero-dimensional ideal. There is a bijection between the set A(7) of its order 
ideals and the set of integral points of the order ideal polytope of I. The bijection is given by 

z eP(7)nZ'^" ^ G{z) := {m e T" | = 1}. 

Proof In fact, we will see that the order ideal polytope is defined exactly to this end. Obviously, the 
order ideals consist of monomials of degree smaller than the dimension d of the factor K[X]/I, so 
restricting to these monomials is permissive. 

First, the integral solutions z of [Condition 3.41 are exactly the 0/1 points, i.e., the characteristic 
vectors of sets of terms 0(z) := {m e T"^_^ | z^ = 1}. 

Second, it is easy to see that [Condition 4.1] means that ^(z) is indeed an order ideal, as whenever 
mi I snd 7772 ^ ^(z), i.e., z^^ = 1, then it follows that z^^ = 1, i.e., e &(z) as well. 

In the third step we will provide an algebraic characterization of Conditions [3.2[ and [3.31 Clearly, 
[Condition 3.2[ can be rewritten to 

mz)\=d. 

We will now show that [Condition 3. 3] is equivalent to 

(^(z))j,n/-{0}, 

i.e., the image of ^(z) is linearly independent in the factor K[X]/I. Thus Conditions [3.2[ and [3.3[ 
together are equivalent to K[X] = I ® (^(z));^, finishing the proof. 
To prove our claim on [Condition 3.3[ we rewrite it similarly as before, 

\Un0(z)\ < dim (UUI)kIi , 

i.e., the size of [/ n ^(z) is at most the dimension of the vector space generated by the image of U in 
the factor K[X]/I. This is obviously necessary for the image of ^(z) to be linearly independent in the 
factor. (Here the size of U does not matter.) 

For sufficiency choose U := 0{z). Then the dimension of the vector space generated by the image of 
0{z) is at least |^(z)|, so the image of 0{z) is independent. □ 

4. Computing general border bases 

We will now turn to actually computing general border bases using Theorem [3. 2[ Clearly, we cannot 
expect to be able to compute a border basis for any order ideal, simply as not every order ideal supports 
such a basis. Having the order ideal polytope at hand something slightly more subtle can be done: We 
can choose the order ideal by optimizing a linear (preference) function c e over the order ideal 
polytope. In his way, we can search for an order ideal with preferred monomials in its support (see 
[Subsection 4.3D . An obvious application is to count the number of border bases that exist for a given 
ideal. 
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Before we can turn to the actual computation though, we will revisit the order ideal polytope and 
provide a reformulation that is suitable for actual computations JSubsection 4.1D . We will then show 
in [Subsection 4.2| how to obtain an ^-border basis for 6 e A(/) where / c K\%\ is a zero-dimensional 
ideal. 

4.1. The order ideal polytope revisited. Classically, the computation of the border basis is performed 
by Gaussian elimination applied to the matrices obtained within LStabSpan. The order ideal for 
which we then compute the border basis is solely determined by the pivoting rule when performing the 
elimination step (i.e., which basis we choose). Term orderings ensure that we obtain an order ideal for 
which a border basis exists (given the correct d). If we would now pivot arbitrarily, which effectively 
means permuting columns, it is not clear that, first, the resulting ideal is an order ideal and, second, 
that it does actually support a border basis. We face the combinatorial problem that when permuting a 
column such that it will end up being an element of 6 we also have to ensure that all its divisors will 
also end up in 6. On the other hand, due to the dimensionality constraint we must choose an exactly 
determined amount of elements that will end up in ^. This problem oi faithful pivoting is captured by 
the order ideal polytope as we will see now. 

Throughout this subsection let M = Ujgf^Mj c K[%] be a finite set of polynomials with degree- 
filtration {Mi I i e N} that generates a zero-dimensional ideal {M)jf[x] Q K[X] such that all p e M; 
have degree i. Furthermore, for each i e N we have an enumeration M; = {pij \ j e [fcj} with fc; e N. 
In the following we assume that L is of the form L = T"^ where d := dimi<C[X]/(M)jf [x]- For the sake 
of clarity, we assume that M is L -stabilized and has in particular a convenient form. One might think 
of M being the output of the LStabSpan procedure, which is then brought into the following reduced 
form: 

Definition 4.1. Let M be a finite set of polynomials of degree at most I for some £ e N. Then M is 
in canonical form if the leading term of any element of M does not occur in the other elements. Here 
we can freely choose leading terms of the polynomials with the only constraint that they have to be 
maximal-degree terms. 

We give a visual interpretation of the definition. The coefficient matrix A & of M is the 

matrix where the rows are the elements of M, the columns are all the monomials of degree at most i, 
and the entries are the coefficients of the terms in the elements of M. We use the convention that for 
terms t^, t2 with deg(ti) > deg(t2) we put column t^ to the left of column Similarly, we put leading 
terms of a polynomial to the left of the other terms. Now M is in canonical form if the matrix A has the 



structure as depicted in Figure 4.1 i.e., it consists of degree blocks and each degree block is maximally 
interreduced. 

The degree blocks correspond to the leading forms of the polynomials in M. Any finite set can be 
brought into canonical form by applying Gaussian elimination and column permutations of terms with 
same degree if necessary. In particular, the output of the LStabSpan procedure can be easily brought 
into this form, which is convenient to work with. 

We are ready to provide a computational reformulation of the definition of the order ideal polytope. 

Lemma 4.2. Let M be T'^^-stabilized and in canonical form and (M)^'^ j {M)^'^~'^ = {t'I^)^^. Let L ■- 



^<d-i assume d = \L\- \M\. Then the order ideal polytope P(M, L) c [0, 1]^ of {M)^[x] is given by 



the system of inequalities in Figure 4.2 



Proof Let I := {M)^^^]- The only difference between the systems in Figure 3.1 and Figure 4.2 is that 
[Condition 3.3l is replaced bv lCondition 4.3[ So we will show their equivalence modulo the other condi- 
tions. 

First, Corollary [2J2] provides = (M)||^j^^ = {M)^. Applying Lemma [iJS]© gives {L)k/{M)k = 
K[%]/I as vector spaces, and so d = \L \ — \M\ = dimK[X]/I. If U' Q L, then via the modular law, we 
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Figure 4.1. canonical form 
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Figure 4.2. Order ideal polytope P(M, L). In lCondition 4.31 the matrix U is the induced 
sub-matrix of M with column monomials only in U . 



have {U'ul)j^ll ^ (U' DM) j^rj {M)y, so we can rewrite [Condition 3.3| to 
(4.4) 2z„<dim([/'uM)j,-|M| 

meL/ 

MU' QL: \U'\ = dim(L)jf - |M|, 

where we have deliberately replaced U with U'. 

We will show that the difference of Conditions |4.2| and |4.3| is equal to [Condition 4.41 with the choice 
U' := L\ U, which has size \U'\ - dim {L)^ - \M\. This will finish the proof. 

Let [/ c L be as specified in [Condition 4.31 and compute the difference of Conditions [4.2[ and [4.3[ We 
obtain 

z^<dim{L)K-\M\-\U\+rKU). 

meL\U 

It is easy to see that \U\ = dim j^lf . We claim that it suffices to show that 

{Mua\u))„ 

rk(L/) = dim — — ; . 
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Indeed, using this we can rewrite the inequality as 

S{L)k 
z„ < dim(L);f - |M| - dim , , , 

meL\U ^ ^ 

{Mua\u))K 

+ dim — — ; 

a \ U)k 

= dim (MU(L\t/))^- |M|, 

which is q4.4D for U' := L\U as claimed. 

We will show now that rk(L/) — dim (M U (I \ U))i^ - dim (L \ U)i^. Let B denote the matrix obtained 
when writing the elements in M as rows and let L index the columns. Clearly, = 0£g^i<Ce£ as a 
vector space and {U')j^ = ©^g^/iCe^ as a sub vector space with U' --LX U. We obtain {L)j^j{U')^ ^ 
®leL\u'^^i- that (M)jf c (L)jf and thus we can consider {M UU') j^j {U') ^ c {L)j^l{U')j^. 

Now dim (^(M U U')kI{u')k) = rk(t/) where f7 is obtained from B by removing the columns in U'. We 
obtain 

rk(LfJ = dim — — = dim — — ; 

and thus the result follows. □ 

4.2. Computing border bases for & e A(7). As the computation of the L-stable span of a set of 
generators M is independent of the actual chosen vector space basis (see Remark |Z61) . we can adapt 
the classical border basis algorithm (cf. [12J) to compute border bases for a general order ideal. We first 
determine the right computational universe L = T" ^ for some d e N such that the associated L-stable 
span M contains all border bases. In a second step we optimize over the order ideal polytope P(M, L) 
and then perform the corresponding basis transformation. We will first introduce the generalized border 
basis algorithm, and then establish its correctness: 

Algorithm 4.3 (Generalized border basis algorithm — BBasis). 
Input: F a finite generating set of a zero-dimensional ideal. 
Output: a border basis of the ideal. 

(1) Let d := max^g^ {deg(/)} and put L := T'^^. 

(2) M = {mi,...,mj := LStabSpan(F,L). 

(3) IfTl^ ^ LT(M) then setd:=d + l and put L := T^^ and go to step ©. 

(4) Set d„i^ := d, d := \L\ - \M\ and L := T^j_^. If d < d„i^ then let M := M-'^. Otherwise let 
M := LStabSpan(F, T-''). 

(5) Choose e A((M)jf [x]) ('^ z ^ P(M-''-\ L) n and = 0{z)). 

(6) Let ^ := BasisTransf orination(M, 0). 

Note that step Q is a convenient way to quickly check whether M is already L -stabilized and suffi- 
ciently large to calculate the dimension K[X]/I. Step Q adjusts d to this dimension and recomputes 
M. 

In this augmented algorithm we added the steps Q, ©, and ©. Step © will be extensively 
discussed in lSubsection 4.3| as there are various ways to determine e A((M)jf[x]) ^^d this is precisely 
one of the main features, i.e., to choose the order ideal more freely. Note that by Theorem 13.21 we 
already know that does support an ^-border basis of (f )j<:[x] (as (F)^ [x] — {^)k[x]') ^nd our task is 
now to actually extract this basis from M. This extraction is performed in step dS]). 

Let A be a matrix representing a set of polynomials where the columns correspond to the monomials 
in some fixed ordering. Let the head (short: Head(a)) of a row a of A be the left-most monomial in the 
matrix representation whose coefficient is non-zero. Note that the notion of head replaces the notion 
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of leading term of a polynomial as we do not (necessarily) have a term ordering anymore. The main 
idea is to reorder the columns of M and then to bring M into a reduced row echelon form such that no 
m e ^ is head of a row of the resulting matrix — a classical basis transformation: 

Lemma 4.4. Let L = TJ^ with £ e N, let M be a finite set of polynomials satisfying {M)^ = (M)^^ [x] n (L)jf 
and let & = {t^, . . . , t^} be an order ideal with d & QL and € A((M)jf[x])- ^^^n Algorithm 14. 51 returns 
the ^-border basis of {M)^^x\- 

Proof First, the algorithm finds I from M. As ^ e A((M)jf [x]) we have that & supports a border basis 
and in particular we have K[X] = (M)^^ [x] ® (^)jf > ^nd hence (I)^^ — {M)i^ ® (^)jf by the modular law. 
We will now show that [Condition (ijj of Definition 12.21 is satisfied. This in turn follows from the fact 



that the algorithm creates every element gj of ^ to have the form gj — bj — ^f^^ OLij t; with a^j e K. 
By construction bj & d and thus 'S is an ^-border basis of (M)jf|-x]. □ 

Algorithm 4.5 (Basis transformation algorithm — BasisTransf ormation). 
Input: M, as in Lemma[ 
Output: as in Lemma [ 

(1) Set £ ■- max^gM deg(m). 

(2) Permute the columns of M such that t is right of m in the matrix representation of M for all 
m e T^^ and t € 0. 

(3) Reduce M using Gaussian elimination: like Algorithm 12.81 but use Head instead of LT and the 
coefficient of the head instead of LC. Let ^' he the result. 

(4) Let ^■-{g^'^' : Head(g) e d 0\. 

(5) Return ^. 

We will show now that Algorithm |4[3] computes an ^-border basis for e A(/). 

Proposition 4.6. Let F = {/i,...,/^} Q K[K] be a finite set of polynomials that generates a zero- 
dimensional ideal I = (F)j<:[x]- Then Algorithm 14.31 computes the -border basis ^ of I for any (chosen) 
e A(/). 

Proof Corollarv 12.121 ensures that upon reaching step Q, we have (M)^^ = I-'^. Via the modular law, 
step sets d to the dimension K[X]/I. The new M and d will also satisfy (M)jf — I-"^. Obviously, 
T"^ contains all order ideals supporting a border basis, i.e., all e A(/) and even the boundary of 
these order ideals. Observe that / = (F)jf[x] = {M)k[x] ^^d thus, by Lemma |4l4l it follows that ^ is 
indeed an ^-border basis of (F)jf[x]- D 

4.3. Computing preferred border bases. Let M and L be as obtained after step @ in Algorithm 
14.31 As shown in Theorem 13.21 and Lemma [4l2l the order ideal polytope P(M-'^~^, L) characterizes all 
admissible order ideals of (M)jf [x] • In particular, this shows that we cannot expect that every order ideal 
supports an ^-border basis of (M)^^ [x] as the characterization is one-to-one. The natural question 
is therefore how to specify which order ideal should be computed, i.e., which z e P(M-''~^, L) n to 
choose. As the coordinates of z are in direct correspondence with the monomials in I we can define: 

Definition 4.7. A preference is a vector c e which assigns a weight to each monomial m ^ L. If 
z € P(M-''~^, L) n Z^, then cz is the score or weight of z. 

As P(M-'^~^, L) c [0, 1]" is a polytope we can compute an order ideal that has maximal score, i.e., 
we can compute Zq e P(M-'^~\ L) n such that 

czo = maxfcz | z e P{M-'^~'^, L) n Z^}, 

with associated order ideal 0(zq). In this sense a preference is an indirect way of specifying an order 
ideal. For certain choices of c e Z^ though it can be hard to compute the order ideal that maximizes 
the score as we will show now. 
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5. Complexity of order ideals 



In this section, we show that finding a weight optimal, admissible order ideal of a zero-dimensional 
ideal given by generators is NP-hard fTheorem I5.2D . The hardness result is unexpected in the sense 
that we merely ask for a nice basis transformation. On the other hand it highlights the crucial role of 
order ideals in describing the combinatorial structure of the ideal. 

From a practical point of view this is not too problematic as, although NP-hard, computing a weight 
optimal order ideal is no harder than actually computing the LStabSpan in general. Further, state- 
of-the-art mixed integer programming solvers that can solve the optimization problem such as scip, 
cplex, or gurobi can handle instance sizes far beyond the point for which the actual border bases can 
be computed (see []4l]). 

We will show the NP-hardness of finding an order ideal that supports a border basis with maximum 
score by a reduction from the |?c-CLiQUE| problem. The latter is well known to be NP-complete (see, e.g., 
fllOII ). Given a graph T = (V, £), recall that a clique C is a subset of V such that for all distinct u, v e C 
we have (u, v) e E. 

fc-CLiQUE: Let F — (V, E) be a graph. Determine whether F contains a clique C of size k. 

We consider the following problem: 

Max-Order ideal: Let M c K[X] be a system of polynomials generating a zero-dimensional ideal and 
let c e be a preference on the monomials. Compute an order ideal &q supporting an ^g'^order 
basis of (M);f [x] with maximum score with respect to c, i.e., 

c(?"H^o)) = max{c(^"H^)) | ^ admissible for {M)k[x]}- 

We will show that for every graph F = (V, £) and k e [\V\] there exists a system of polynomials 
P\v\,k — K[Xy I V e V] spanning a zero-dimensional ideal such that solving the lMAX-ORDER iDEALl prob- 
lem for solves the |fc-CLiQUE| problem for F. For this, we construct an ideal encoding all /c-cliques 
of the complete graph on n vertices: Let n e N and k e [n] and define 

Pn,k ■■= {Vj lie [n-?c]}UT^3 

with Vj := Xiie[n] ^^-"-i- Consider the ideal generated by F„ ^.. We show that its order ideals are in 
one-to-one correspondence with the fc-element subsets of the set of n variables x^, . . . , as stated in 
the following lemma. 

Lemma 5.1. Let n e N and k e [n]. Then F^j^ generates a zero-dimensional ideal such that 6 e 

^((^".'^)jf[x]) ^"^^ °^^y - ^=1 i^^t^ l^'^l = ^' ^ I ^'y ^ ^"H, and 6=^ = 

for all l>3. 

Proof We start by providing an explicit representation of the factor ring iC[xi, . . . , x^]^ (^".'=)k[x]' 
The polynomials Vj with i e [n — k] are homogeneous of degree one. The coefficient matrix A := 
(Vj)jg|-„_j,] is actually a Vandermonde matrix and in particular every square submatrix of A is invertible. 
In particular, removing the columns belonging to any k variables of {x^, . . . , x„} from A, without loss of 
generality say, x^, . . . , Xj,, the resulting square submatrix is invertible. Thus, {x^, . . . , Xj,, v^, . . . , v^.j.} is 
a basis for the homogeneous polynomials of degree one, and hence the ring is also a polynomial ring 
in these variables. In particular, 

{^)k = K[x„ . . . , X j/ = K[x„ . . . , X,]/ (Ti,)^^^^ . 

As substitution preserves degrees, homogeneity, etc., it follows, that any order ideal has to have 1 \ = 
k, = {xy I x,y e 0='^}, and = for all £ > 3. 
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The other direction follows immediately as each order ideal 6 with |^"^| = fc, 0^^ — {xy | x,y € 
0"^}, and = for alH > 3 is actually an order ideal such that = {Pn,k) j^^j^-^ ® (^)if' 

as above. □ 

Note that the order ideals of F^j^ indeed correspond to the fc-cliques of the complete graph on n 
vertices: If ^ e A(F„jt), then — and Xj Xj^ e if and only if Xj^, e 0"^. If we 

now remove all elements of the form x,^ with x, e 0^^, and there are k of those, then 

fc(fc-l) 

the size of a fc-clique. It is worthwhile to observe that all admissible order ideals supporting a border 
basis of {Pn,k) degree-compatible as f„ j, is a homogeneous system. Thus the reduction pre- 
sented in [4] would also establish NP-hardness in the general case. We will now present a more direct 
reduction not relying on the degree-compatible case. 

Theorem 5.2. IMax-Order idealU s NP-hard. 



Proof. The proof is by a reduction from the NP-hard [A:-Clique[ Let the graph F = (V,£) with n := \V\ 
and ?c e [n] be an instance of lfc-CLiQui] We consider M := F^^. and define c e Z^i^ via 



c = 



1, if m = x^Xy and either (u, v) e £ or u = v; 
0, otherwise. 



By Lemma [57T1 the order ideals of (M)jf[x] are in one-to-one correspondence with the fc-cliques of the 
complete graph on n vertices. A vertex v e V is represented by the degree-one monomial x^ and each 
edge (li, v) e £ is represented by the degree-two monomial x^x^. Thus, the score of an order ideal is 
the sum of k and the number of edges in F between the k vertices of the clique. 

Clearly, there exists an admissible order ideal with score MMil if and only if F contains a clique 
of size k. We obtain that IMax-Order idealI solves |fc-CLiQUE| and so the former has to be NP-hard. □ 
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